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We provide provide a detailed study of biasless coherent transport of singlet electron pairs in 
one-dimensional (ID) channels induced by electron-electron interactions that are time-varying in 
certain spatially localized regions of the channel. When the time variation is cyclic, the mechanism 
is analogous to the adiabatic quantum pumping of charge and spin previously studied. However, 
the presence of interactions that vary only in localized regions of space requires an intrinsically 
two-body description which is irreducible to the ID single particle scattering matrix elements that 
are sufficient to describe quantum pumping of charge and spin. Here we derive a generalized theory 
for the pumping of such interacting pairs starting from first principles. We show that the standard 
description of charge pumping is contained within our more broadly applicable expressions. We 
then apply our general results to a concrete lattice model and obtain an exact analytical expression 
for the pumped singlet current. We further demonstrate that such a model can be implemented 
with a chain of currently available quantum dots with certain minor modifications that we suggest; 
we present a detailed numerical feasibility analysis of the characteristics of such experimentally 
realizable quantum dots, showing that the requirements for a measurable pumped singlet current 
are within experimental range. 



PACS numbers: 72.25.Pn, 03.67.Mn, 73.23.-b,73.63.-b 



I. INTRODUCTION 



Pumping mechanisms have attracted research at- 
tention in recent years, both experimental h&SdaO&iLIL 
and theoretical 9 ! 10 ! 11 ! 12 ^ 3 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 ! 23 , be- 
cause they provide a remarkable alternative means of 
generating current in nanostructures. Such mechanisms 
use time varying parameters or potentials to achieve the 
flow, rather than the application of a bias. While the 
basic idea behind such a pump follows from classical in- 
tuition, in nanostructures quantum mechanics can in- 
troduce novel features that make "quantum pumping" 
a rich and non-trivial process. For example, a quan- 
tum dot with time varying couplings to pinched off leads 
will pump electrons in a process that is analogous to the 
classical pumping of water through a lock in a canal^i 2 -; 
here Coulomb blockade plays the role of gravity, limit- 
ing the amount of charge that can flow through the dot 
each pumping cycle. However, in an open quantum dot 
system, time varying parameters lead to changes in the 
electron wavefunctio n 24 ! 25 , leading to modulation of the 
spatial probability density that results in a net transfer 
of particle density through the region of time variation. 
This quantum mechanical process is subtle and cannot 
be understood with classical intuition alone. 

The concept of quantum pumping is a very broad one, 
generically relevant when a channel has time varying pa- 
rameters. As a result, one can imagine using pumping 
to transport any quantity associated with quantum me- 
chanical wavefunctions that can be varied in time with 
relative ease. Charge 6 -* 7 -! 9 - ! 12 ' 15 ! 16 and spirt 5 - ! 11 ' 14 ! 20 have 
been primarily considered for pumping. These quantities 
can be associated with single particles. Their pumping is 
most often computed within an independent particle ap- 
proximation, although there have been studies that have 



computed their pumping in a strongly correlated system 
like a Luttinger liqui d 26 ' 27 . 

In a recent paper—, we proposed the quantum pump- 
ing of a two particle quantity - electron singlet pairs. 
This is relevant from an applications standpoint as well 
as a fundamental physics point of view. A pair of elec- 
trons in a singlet state is spin-entangled 2 ^, and coher- 
ent transport of entangled states is essential for quantum 
information 29 . Describing the pumping of a two parti- 
cle quantity is quite different from the usual quantum 
pumping theory of single particle quantities. As one of 
the two main topics of this paper, we derive a theory that 
is capable of describing singlet pumping and also incor- 
porates standard quantum pumping cases. The inher- 
ent presence of two body interaction in pumping singlet 
pairs makes the theoretical description significantly more 
involved than the commonly used models based upon sin- 
gle particle scattering matrix elements 9 -. Thus a primary 
goal of the paper is to provide a detailed analysis of the 
necessary theory for a singlet pump. 

The second topic of this paper is a study of the experi- 
mental feasibility of a quantum singlet pump based upon 
a detailed numerical analysis of an experimentally realiz- 
able system involving quantum dots. We find that, with 
laboratory systems already available, the implementation 
of a singlet pump could be possible. 

The paper has roughly three parts. The first part in- 
volves the derivation of our approach to quantum pump- 
ing and comprises of Sec. II, III, and IV. Section II de- 
fines a singlet current and interacting two particle states. 
Section III develops the theory of a singlet pump starting 
from an adiabatic perturbation expansion. These sec- 
tions rely on Green's functions theorems appearing in 
Appendices A and B. Section IV demonstrates that our 
results recover the established theory when one considers 
quantum charge pumping of single electrons. The second 
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part of the paper appears in Section V which describes 
a physical model based upon a two particle analog of a 
turnstile, involving a chain of quantum dots. An ana- 
lytic expression for the singlet current in that model is 
derived using the Green's function theorems in the ap- 
pendices once again. The third and last part compris- 
ing of Sec. VI, VII and VIII is devoted to a feasibility 
analysis of a proposed experimental implementation of a 
singlet pump based upon the physical model presented 
in Sec. V. Section VI suggests use of a specific design of 
quantum dot similar to an experimentally realized dot, 
and the spatial potential energy profile of an electron in 
the dot is numerically computed. Section VII evaluates 
the energy of two electrons in such a dot, and finally 
Sec. VIII demonstrates that such a dot has the features 
necessary for generating a pure singlet current while sup- 
pressing current of single particles and triplet pairs. 



II. SINGLET CURRENT 

A. Definition of the current 

The quantum mechanical current density is generally 
defined from the continuity equation by considering the 
time variation of the single particle density. In order to 
discuss quantum pumping of singlets, we need to have an 
appropriate definition of a singlet current. We thus begin 
with another continuity equation. The probability that, 
within a one-dimensional system, there is one electron 
at position and spin X\ = (xi,ai) and another at X 2 is 

^(X u t)^(X 2 ,t) 1 p(X 2 ,t)^(X 1 ,t)y Here, ft(X u t) = 
e iHt/h^\^ Xi y r iHt/h and ^t(Xi) creates an electron at 
X\ . Here and henceforth, all expectation values are taken 
with respect to the N particle state of the system \N). 
Averaging over the position x 2 of the second electron, we 
find from the Schrodingcr equation 



d_ 

dt 



dx 2 U\X 1 ,t)^(X 2 ,t)^(X 2 ,t)^(X 1 ,t) 



d 
dxi 



h I ( d d 
2mi J 2 \dxi dx[ 



(1) 



x ( ft(X u t)ft(X 2 , t)i>(X 2 , t)j>{X[,t) 



X 1= X[ 



We are led to the definition of the spin-specific two- 
particle current density 

J(xi,ai,a 2 ,t) (2) 
eh f 

= — / dx 2 Im{d Xl p 2 (X 1 ,X 2 ;X[,X 2 ,t)} Xi=x , i 



in terms of the two particle reduced density matrix 



Pz{Xi, X 2 ; X[, X 2 ,t) — 

(^\X 1 ,t)^(X 2 ,t)i>(X^,t)i > (X[,t)). (3) 
A summation over the two spins would yield a current 
analogous to the usual definition of current but for the 
flow of pairs of particles rather than single particles; for a 
single pair of particles that are in momentum eigenstatcs 
we would obtain the average current density of the pair. 

We can expand the the two particle density matrix in 
terms of energy eigenstates discussed in Appendix C, 



p^X^X^X'^X^t) (4) 
= J dEf{E)^ E {X u X 2 ,t)^l E {X[,X' 2 ,t) 



where E signifies the available energy of a pair of particles 
and f(E) is a distribution function for the pair, which we 
take to depend only on the total energy E of the pair. As 
we demonstrate in Sec. III.C, this is consistent with the 
fundamental physical requirement that no current flows 
in the absence of pumping or bias. 

In the absence of interaction, the two particle density 
matrix can be separated into singlet and triplets spin 
subspaces. The pair interactions that we consider do not 
affect spin by assumption, hence such a separation would 
continue to apply. Furthermore, the energy eigenstates 
would also be spin eignstates, so that we factorize the 
states ^E(Xi,X 2 ,t) = ^ E (xi,x 2 ,t)xu, where \v is a 
singlet xs or a triplet Xt„= ,±i s P m state, and effects of 
the time-dependent potential is felt only by the spatial 
part (%i , x 2 , t) . Thus we write the two particle density 
matrix as 



P2 = Ps ' 



XsXg "I" 



cr=0,±l 



(5) 



where p s = ^* EtS (xi,x 2 ,t)^ E ,s{xi,x 2 ,t) and p Ta = 
ty* E T ^(x\, x 2l t)^E.T a (xi, x 2 , t) denote the spatial com- 
ponents. Energy dependence on the spin states is im- 
plicit. We can then define the current density for each 
spin subspace. For singlets we have 
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eh f 

Js(x,t) = — / dx 2 lm{d xi p s {x 1 ,x 2 ;x' 1 x 2 ,t)} , (6) 
m J 1 1 

= & ±J dEf(E) jdx 2 j d ±j ^ ^^ + ^-E)l m {^^ 1 , fc3 (x 1) , 2!Wl , fe2 (x' l! x 2 ,t)} si=x , 



r 



where we have parameterized the spatially symmetric sin- 
glet states ^E,s(xi,x 2 ,t) = $ kl ,k2{xi,x 2 ,t) by single- 
particle momenta k\ and k 2 corresponding to the nonin- 
teracting momentum eigenstates from which these states 
can be generated as we show in the following subsec- 
tion. The symbol f(E) denotes the state occupation only 
within the singlet subspace. 

We note that Js gives the flow of probability of finding 
one member of a singlet at X\ irrespective of the location 
of the other member. If we average over the position of 
the other particle, we will obtain the expectation of the 
current over the length of the system. 



B. Interacting two particle states 

Suppose our system has the time-independent many- 
body Hamiltonian, 

H = J dX^(X)h(x)i>(X) (7) 
+ J dX J dX'^(X)^(X')V(x,x')^{X')^{X) 

where h(x) — —^-^i + W(x) is a first quantized single 
particle Hamiltonian with some external potential W{x). 
The equation of motion for ^ kl>k2 (xx,X2,t) can be de- 
rived directly from the Schrodinger equation. One finds 
that with suitable approximations (see Appendix C) a 
two-particle equation arises 

d 

~ h ( x V ~ M 3 ^) ~ V(xi,x 2 )]^ klM (x ll x 2 ,t) = 

at 

[Eih,^) - H Q - V(x 1 ,x 2 )]^ klM (x 1 ,x 2 ,t) = (8) 

where Ho = h(xi) + h(x 2 ) and ^k 1 ,k 2 (xi,x 2 ,t) 
has the trivial time dependence ^ kl ,k 2 i x i, x 2 , t) = 
e- iE ^ k ^ h ^ kuk2 ( Xl ,x 2 ,t). 

If the one body potential and two body interaction 
were absent, the solution would simply take the form of 
free singlet states 

®k(x) = -^[0 fcl (xi)0 fc2 (x 2 ) + <t> kl { x 2)4>k 2 {xij\ (9) 

where 4>k denotes a single particle momentum eigenstate 
with momentum hk, and where we have introduced the 
notation x = {xi,x 2 } and k = {ki,k 2 }- Even in the 
presence of a one-body potential, this form still holds 



except that the single particle states would be defined by 
h(x)(j>k{x) = E k <fik(x). 

In the presence of the two body interaction V(x) — 
V{xi,x 2 ), the interacting or scattering singlet states ^ 
can be expressed in terms of the free state $ by the 
Lippmann-Schwinger (LS) equation 

* s (i) = + j dx'G(x,x'-,E)V(x')$- k (x')(10) 

where the full retarded two particle Green's function sat- 
isfies 

[E-Hq- V(x)]G(x, x'; E) = 5{x - x'). 

A few comments about expression (|10[) are in order. 
First, due to the two-body interaction, the scattering 
state generally involves a range of momenta, so the sub- 
script k there simply serves as a book-keeping label 
to indicate the non- interacting state that generates it. 
Secondly, although the wavefunctions are singlet func- 
tions, the expression involves the unsymmetrized full re- 
tarded two particle Green's function G(x, x'; E). Finally, 
Eq (fTO)) applies to any two particle wavefunction, regard- 
less of its symmetry. If we choose a symmetric, spin- 
independent potential, the symmetry of the free states is 
preserved by the scattering and & k (x) will have whatever 
symmetry that ^ k (x) has. In our case, we are interested 
in the evolution of a spatially symmetric singlet state, 
so we apply (fTU|) to this case. It is not as interesting to 
consider the pumping of triplet wavefunctions since the 
current could consist of an unspecified superposition of 
all three different triplet spin states. 

III. PUMPED SINGLET CURRENT 

A. Pumping via Interaction 

Since the equation of motion (fT0)l does not inherently 
distinguish between singlet and triplet states, if we wish 
to pump only singlets we need to find a potential that 
only affects singlets. Obviously, manipulating the ex- 
ternal potential W(x) in ([7]) would affect single particle 
states as well as singlets. Thus, we focus on the two- 
particle interaction V(x, x'). Within our one-dimensional 
system we consider a time-varying two particle interac- 
tion that exists at two specific localized regions, as shown 
in Fig. [T] Two electrons interact only when they are both 
in one such region. If either electron is outside the re- 
gions, there is no interaction, and we assume that the 
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interaction vanishes between an electron at one region 
and an electron at the other region. Physically, if each 
region is sufficiently well localized in space, and if the 
one-dimensional system consists of a single channel, the 
Pauli Exclusion principle will disallow two electrons in a 
triplet state from both occupying the same region. As 
a result, electrons in a triplet will not feel the interac- 
tion. Only singlets will feel its effects. The specifics of 
how such a localized interaction could be implemented 
are discussed in Sec. V. For our derivations below, we 
will assume that it is exactly true that the interaction 
V(x , t) at each site affects only singlets, and that triplets 
are completely unaffected by it. 



B. Adiabatic Perturbation 

The most commonly used description of quantum 
pumping employs the well-known Brouwer formula- that 
derives directly from the Landauer-Buttiker formalism. 
Another common description uses Flouquet theor y 18 ' 19 . 
Since those approaches rely upon single particle scatter- 
ing matrix in ID, not appropriate for describing inter- 
acting pairs, we have taken an a priori approach rooted 
in adiabatic perturbation as used by Thouless^ in his 
original paper on quantum pumping. That approach has 
been used in some recent papers^ as well. This will 
merely constitute our point of departure, because we have 
the significant complication of interacting particles with 
pumping driven by the interaction, while the previous 
studies were derived in a single-particle picture. 

The singlet current in Eq. ([6]) is defined in terms of a 
time-dependent singlet two particle state &k(x). There- 
fore the determination of the pumped current reduces to 
describing the evolution of this singlet state in the pres- 
ence of a time-varying interaction V(x,t) 



H = H + V(x,t) 



(11) 



in Eq. ([8]). Here, the free Hamiltonian Ho is still time- 
independent and its two-body energy eigenstates are still 
of the form shown in Eq. @ comprised of symmetric 
combinations of single particle eigenstates. 

For an adiabatic process, the potential, or the two- 
particle interaction in our case, is assumed to vary slowly 



(a) 

Contact 



V(xi, x 2 ; t) 



Contact 



FIG. 1: (a) Schematic figure of a one dimensional system in 
which the two-body interaction V(xi, X2 \ t) acts when x\ and 
X2 are both within a finite interval near —I or when both are 
within a finite interval near I. 



compared to the time spent by the particles in the re- 
gion of the potential, i.e. the dwell time is much shorter 
than the time for variations in the potentiate. As a 
result, it is appropriate to insert the ansatz ^^(x,t) = 
e -ijdtE(k,t)/n^,^,t) into Eq. ©, we are left with 

*) = [Ho + V(x) - E(k, t). (12) 

Henceforth we will drop the tilde symbol on the phase 
e -*/ dtE(k,t)/h ^hat distinguishes \& from '5 does not af- 
fect the current ^ in any case. To find a zcroth order 
solution, we neglect the time derivative and solve the in- 
stantaneous equation 



[E0t,t)-H o -V(x,t)]*\{x) = 



(13) 



where the time t is simply a parameter. This equation 
on inversion becomes the Lippman-Schwinger equation, 
Eq (flU]) , for the instantaneous scattering state To 
find corrections, we evaluate the instantaneous Green's 
function 

[E-H -V(x,t)]G t (x,x';E) = 5{x-x), (14) 
and note that the exact solution to (11211 is 



* s (2,t) = *|(S)- « J tte'G t (i > x / ;£?)^*fc(2 / ,t).(15) 

Iterating this equation, we can compute corrections to 
the zeroth order solution ty%(x). To first order, 

= *|(.x) + A%(^) (16) 

Higher orders are left out here because of the assumption 
of adiabaticity. Taking the derivatives of Eqs (|13I14[) with 
respect to time yields the following important identities 
for the time derivatives of the instantaneous functions 



= / dx'G t (x,x';E)V(x',t)$>l ; (x') (17) 
G*(x,x';E) = / dx"G i (x,x , ;E)V(x",t)G t (x",x' ; E) 



which allow us to write the second term in Eq. (|16|h the 
term first order in the time derivative) as 



A%(sc,t) = -ih J M'G t {x,x';E)^! t (x') (18) 
= -in f dx' f dx"G* (x, x'\ E)G t (x, x"; E)V(x")^ t (x") 



while the zeroeth order term is simply (x) determined 
by the Lippmann-Schwinger equation (| lOf) . 

Henceforth, we sometimes simplify notation by writing 
G(x,x') instead of G (x, x'\ E). 
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C. Zeroth order: No Spontaneous Current 



At the zeroth order in the time-dependence, at each 
instant the system is unaware of the fact that the po- 
tential is changing. Since there is no bias either, the 
current should vanish, because a non-vanishing current 
at this order would essentially be a spontaneous current, 
which is unphysical. In order to demonstrate that our 
definition of the singlet current passes this crucial test, 
we explicitly evaluate the current at the zeroth order in 
the time-derivative, which is given by 



Jo 



Jo = 



*«(^+^- B ) Im <**" (i)s *'*' (i) } 



Using the following identity for the free two particle 
Green's function 

Im{G (5, x'; E)} =-J. dk5 - e) $fc(x)$* (x'\20) 



together with the Lippman-Schwinger Eq. 
reduce the expression (fll?|) to 



dE f(E) / dx' / dx"V(x')V{x") (22) 



eh 
2irm 

x Im|iG(s";x') J dx 2 [G {x' , x)d xl G (x; x")] 
- iG*(x"; x') [ dx 2 [G* (x; x")d Xl G* (x, x')] 



The kernel of the integral is of the form F x 
Igg(x%,x',x") evaluated in Eq. (jBlip in Appendix B 
(here F = iG(x', x"; E)). Such an integral was shown 
to vanish for |xi| > \x[ |, \x"\. This applies here since the 
current is measured outside the region of interaction and 
the variables x[ and x'{ are associated with the localized 
interaction factors V(x') and V(x"). Therefore the net 
current at the zeroeth order vanishes, Jq = . 

D. First order: Adiabatic Pumped Current 

The pumped current to first adiabatic order is a bilin- 
ear form involving the zeroth order wavefunction ^I(S) 
and the first order wavefunction correction A^/j,(a;,i) 
fromEq. (JSJ|: 



Jo 



-eh 



dEf{E) / dx 2 lm[d Xl lm{G (x,x)} 



-h / dx" {d Xl G(x;x")}V(x")Jm{G {x" ,x)} 

dx' G*(x;x')V(x')d Xl Im{G (x,x')} 
- I dx' [ dx" {G*(x;x')V(x')} 



d Xl {G(x;x")V(x")Im{G (x",x')}] (21) 

The first term vanishes immediately because it involves 
the imaginary part of a real number. Then after some 
manipulations where we use the Lippmann-Schwinger 
equation, the Dyson equation G = Go + GoVG, along 
with the reciprocity property of the Green's functions, 
the expression for the current at the zeroeth order can 
be reduced to 



xIm{^f(x)d Xl A^^x,t)+A^l(x,t)d Xl ^l(x)}(23) 

Using the identity in Eq. (|20|) and by multiple use of 
the Lippmann-Schwinger equation and the Dyson equa- 
tion, in a straightforward but lengthy calculation, the 
above expression can be reduced to 



Ji = 



27TTO 



dEf(E) / dx' / dx"V{x") / dx 2 (24) 



x \{G*{x",x) - G(x",x)}d Xl {G{x,x')G(x'x")} 
{G* (x, x')G* {x'x")} d Xl {G*(x,x") - G{x,x")}\ . 

The identity / dx'G(x, x')G(x' , x") = -d E G(x,x"), 
(Appendix A) then transforms this to 



Ji = 

•ha = 
Jlh = 



Jla + Jlfc 

2irm I oE 



dx"V(x") / dx 2 lm{G*(x",x;E)d xl G(x,x";E)} 



2itvn 



dE f(E) / dx' / dx"V(x")Im{G*(x',x";E) / dx 2 G*(x,x';E)d Xl G*(x,x";E) 



+G(x',x";E) I dx 2 G(x",x;E)d Xl G(x,x';E)} 
I 



(25) 



We will now show that for any observation point x\ greater than the region of non-vanishing U, the term J\ 
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would not contribute in a periodically varying pumping 
cycle. We should first point out that the pair of terms in 
Ju are not complex conjugates, due to the assymmetry 
with respect to the exchange of x' <-> x" , so we cannot 



J 



argue that taking the imaginary part makes it vanish. 
Using the Dyson equation we can make the following ex- 
pansion 



lm{G*(x', x"; E) J dx 2 G* (x, x'; E)d xl G*(x, x";E) + G(x', x"; E) J dx 2 G{x" ', x; E)d xl G(x, x'\ E)} 

lm{G*(x', x")fdx 2 G*(x, x')d Xl G*(x, x") + G(x', x")Jdx 2 G (x, x")d Xl G (x, x')} 
J dy Im {G*(x\ x")V(y)G* % x") J dx 2 G* Q (x, x')d xi G* (x, y) + G(x', x")V(y)G(y, x") J dx 2 G (x, y)d Xl G (x, x')} 
Jdz Im{G*[x\x' / )G*(z,x')V(z)Jdx 2 G^[x,z)d Xl G^(x,x'') + G{x\^ 
Jdzjdy Im {G* (x', x"){V(y)G* (y, x')V(z)G*(z, x")} J dx 2 G* (x, y)d Xl G* (x, z) 

+G(x',x"){U(y)G(y,x")U(z)G(z,i')}/^2G (x,y)a :!;i Go(S,z)} (26) 

I 



All four terms on the right hand side of the expression 
above involve integrals of the form F x Igg( x i> & , m 
Eq. (|B1 1[) in Appendix B, the function F being different 
in each term. Using the properties of that integral in the 
Appendix we see easily that the last two terms would 
vanish for x\ outside the interaction region. In fact if 
the observation point |xi| — > oo, which would guarantee 
|ari| > \xi\ , by the same argument the first two terms 
would be zero as well. However since the variable x[ is 
not associated with any interaction factor V, in principle 
it is not bounded and we have to allow for \xi\ < \x[\; 
in that case Eq (|B12p implies that for observation points 
outside the interaction region \xi \ > \x'{\ , we can reduce 
the first two terms on the right hand side (which provide 
the only non-vanishing contributions) to 



2-^{6(x[ - xx)6(xi - x'{) - 6(x'l - x 1 )9{x 1 - x[)} 
xIm[G(x', x"){G (x', x") + JdyG (x, y)V(y)G(y, x")}} 



= 2^lm[G(x',x")G(x',x")] 
x{e(x[ - xi)6(xx) - 6(-x 1 )d(x 1 - x[)} 



(27) 



Inserting this into the second line of the expression for 
Jib in Eq |25|) and using the identity in Eq (|17|) we get 



the following 

1 Q r , , , 
Jib = -tt / dx {6{x x -xi)6{xi) - 6(-xi)d(xi -%)} 
7T Ot J 

x J dE f{E)lm{G{x',x';E)} (28) 

= -^-Jdx'p(x') 

x{6{x[ - xi)9{xi) - 9(-Xi)d(x x - x[)} 

In the last step we use the definition of the instantaneous 
two-dimensional density of states iIm{G(i', x 1 ; E)} = 
—p(x', E) and did the integral over energy to get the lo- 
cal density p(x'). As we see, this is a total time derivative 
that would vanish on integrating over a full period. In 
any case for an observation point \x'\ 3> even without 
doing the time-integral, this term can be made negligi- 
ble, and actually in the single particle picture used to 
describe charge pumping, the equivalent of this term is 
identically zero since the observation point is taken to be 
asymptotically far from the region of interaction. Thus 
we can now write a compact expression for the net con- 
tribution to the adiabatic pumped singlet current to first 
order for periodically varying two particle interaction, 



Ji(t) 



eh 2 
277m 



dE f(E) — 
M 'dE 



dx'V(x') / dx 2 lm{G*(x',x;E)d Xl G{x,x';E)} 



(29) 



r 



IV. VALIDITY FOR THE NON-INTERACTING 
CASE: CHARGE PUMPING 

Our result has the advantage that it also describes 
adiabatic quantum pumping where the independent par- 



ticle description is assumed. We simply interpret the 
functions as single-particle objects, i.e. the green's func- 
tions are single particle green functions, G(x' , x"; E) — > 
g(x' , x" ; E), and we then do not have the integral over 



7 



X2, so that 



(30) 



dx'V{x')\m{g*{x', x; E)d x g{x, x'; E)} 



Noting that the single particle Green's function has the 
asymptotic form 



g(x,x";E) 



ikfi 



■e tkx ^*{x') 



(31) 



where k = \/2mE/h , one obtains 



J(x,t) 



2nh 



3 , dEf(E)^ilw\V\i>)}(:V2) 



which agrees with the expression in Ref. [16j on noting 
that in that paper the states are normalized by Vk. That 
expression in turn has been shown to be equivalent to the 
Brouwer formula 9 . 



V. PHYSICAL MODEL FOR A SINGLET 
PUMP: A LATTICE OF QUANTUM DOTS 

In theoretical studies of quantum pumping, a turnstile 
modell^ has been used, wherein delta function potentials 
exist at two points in a ID system. The cyclic variation 
of the strength of the potentials serves as the pumping 
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FIG. 2: (a) A discrete model for a singlet pump comprising 
of a chain of quantum dots where the interaction between 
a singlet pair is non-vanishing only when both particles are 
together at one of the two shaded dots, (b) The continuum 
version of the same model where the pair interaction is non- 
zero only at the two points x = ±a indicated by spikes. The 
strengths of the two-particle interaction at the two sites are 
U±(t) which are the time- varying parameters. At any instant 
the two parameters can be of different magnitude indicated 
by different shades in the discrete model and different heights 
of the spikes in the continuum model (the short dotted line is 
simply a marker for the origin x — 0). 



cycle. The region between the two spikes can be inter- 
preted as a scattering region, such as a quantum dot, and 
the strength of the potentials a measure of the barriers 
segregating the scattering region from the rest of the sys- 
tem (e.g. the leads). We can adapt the turnstile model, 
proposing a singlet pump in which delta function external 
potentials acting at the two points are replaced with the 
electron- electron interaction acting only at those points. 
In a continuum of positions, this model of extremely lo- 
calized interactions seems somewhat unrealistic. How- 
ever, in a ID tight-binding lattice, we would have inter- 
actions at two lattice sites, which seems physically rea- 
sonable and appropriate for the application of our results. 
Physically such a lattice model could be implemented, for 
instance, by a chain of quantum dots where precise volt- 
age controls produce only time-varying interaction be- 
tween the members of an electron pair when both parti- 
cles occupy one of two specific dots. A schematic of the 
model is shown in Fig. [51 

Formally the expressions in the discrete case are identi- 
cal with those of the continuum case, as we have verified 
explicitly-^: The Green's functions in the continuum case 
need only be reinterpreted as discrete Green's functions 
and spatial integrals need only be replaced by spatial 
sums. Since the general derivation in this paper has been 
in a continuum model, we will derive the results for the 
turnstile model also in a continuum form. 

The interaction is 



V{x,t) = U-{t)6(x + a) + U+(t)6(x - a) 



(33) 



where ±a = ±{a, a}. The strength of the delta func- 
tions U- and U+ are the two time-dependent pumping 
parameters. The instantaneous pumped singlet current 
in Eq (|2"5)) then becomes 



eh 2 
2-nra 

XC?£;Im 



dEf(E) ]T U±{t) 

y=±a 



dx 2 G*(x,y)d Xl G(x,y) 



(34) 



The full instantaneous Green's function for this system 
can be written compactly as 

G(x,y) = G (x,y) + F-(y)G Q (x,-a) + F+(y)G (x,a) 
with the coefficients being 

"T ± G (±a,y) +T ± T T G a (2a)G Q (Ta,y) 



F±(V) = 



l-T T T±G {2a)G (2a) 



(35) 



We have introduced the T-matrix for a single Hubbard 
interaction (a single term in Eq. (|33p b 



T ± (t) = [U± x (t) + G (0)]- 1 . 



(36) 



We have also introduced a shorthand for the free Green's 
functions: Since the free Green's functions depend only 
on the absolute difference of the coordinates, Go(x,y) = 
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Go(\x — y\), we define Go(2a) = Go(a, — a) = Go(— a, a) 
and Go(0) = Go(a,a) — Go(— a, — a). Using the expan- 
sion Eq. (|35f we can evaluate the coordinate integrals in 
Eq §M§ to be 

— / dx 2 G*(x,±a)— G(x,±a) = (37) 
J a Xl 

-Im{G o (0)} [1 + [F ± (a) + Fl(a)} + \F T (a)\ 2 + \F T (a)\ 2 ] 

-lm{F T (a)G 0R (2a) - F*(S)G* 0R (2a)} 

-lm{FUS)F T (d)G R(2a) - F*(a)F±(a)G* 0R (2a)} 

In order to obtain this form we used the expansion 



Eq (|35|) above to reduce the above integral to a sum 
of integrals of the form shown in Eqs. (jBlft and (|B4|) 
considered in the Appendix B. Then we used their re- 
spected evaluated forms given in Eqs. (|B5j) and (|B6[) . 
What is noteworthy is that the expressions depend on 
Gor(2cl), the regular part of the two-particle Green's 
function which has no singularity. This is discussed in the 
Appendix in Eq. (|A6|) After carrying through a lengthy 
but straightforward set of algebraic manipulations, of the 
above expression, we insert the result into Eq (|34|) to get 
the following expression for the singlets pumped in a sin- 
gle cycle of period r 



J 



s(r) = ^ jdtjdEf{E)—Y.U±{t) 



dE 



\T ± (t)\ 2 [lm{G o (0)}(l + \T T (t)G (2a)\ 2 ) ± 2hn{T T (t)G (2a)G± R (2a)}] 
U±(t)i\l - T T (t)T±(t)G (2d)G (2a)\2 



(38) 



where Gq R (2o) 



G 0R {2a) and G 0R (2a) 



G* R (2a). 



Thus in order to evaluate the current we only need to 
evaluate the free two-particle lattice Green's function, 
for only two specific arguments Gq(Q) = Go (a, a; E) and 
Go(2a) = Go(a, —a; E). Exact analytical forms exist 
for the lattice Green's functions, Go, in terms of ellip- 
tic integrals^. 

As mentioned above the expression for the discrete 
lattice is identical to this with the replacement of the 
Green's functions by their discrete equivalents Go(0) — > 
Go(fh,fh; E) and Go(2a) — > Go(fh, —fh; E) where ±m 
correspond to the lattice sites where the time varying in- 
teraction occurs. In Fig. [2](a) , we take m = 1. The effect 
of changing the value of m on the number singlet pairs 
pumped per cycle was presented in an earlier paper—. 



VI. CHARACTERISTICS OF A QUANTUM 
DOT FOR TIME- VARYING INTERACTION 

While the physical model comprising of a chain of 
quantum dots as presented above yields a measurable 
singlet current in theory, the actual experimental imple- 
mentation requires a feasibility study involving realistic 
quantum dots. We will now suggest a specific configu- 
ration for an experimentally realizable quantum dot to 
tailor it for singlet pumping. In the rest of the paper, 
we will establish that the characteristics of our proposed 
configuration will allow the variation and manipulations 
to satisfy the requirements of a singlet pump. 

To maintain close contact with experiments, 
we specifically consider a cylindrically symmetric 
GaAs/Al r Ga-|_. T A s q uantum dot used in experiments 
reported in Ref. 132- To adapt the dot for singlet 
pumping, we propose a single modification: instead of 
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FIG. 3: (a) Schematic figure of the quantum dot used in Refs. 
[3^ ] and [3^. Here V is the potential energy at the top gate, 
(b) Our proposed modification of that structure with two con- 
centric metal gates, with independently controllable voltages; 
Vin is the net potential energy of an electron at the top edge 
of the blocking barrier (Vi n = — e x where $i n is the 
voltage) right below the inner gate (darker shade). Likewise 
Vout is the potential energy of an electron at the base of the 
outer gate (lighter shade). 



being one single disk, the top gate needs to be made 
of an inner disk of radius Ri n and a concentric outer 
annulus with inner and outer radii Ri n and R ou t, so 
that the voltages Vi n and V ou t on the two pieces can be 
controlled independently. This provides two independent 
parameters for manipulating the shape and depth of the 
lateral potential. A schematic cut-away view is shown 
in Fig. O in the original experiment (a) and with our 
proposed modification (b). 

Along the direction of the z-axis (Fig. [3]), electron con- 
finement is provided by the difference in conduction band 
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edge energies between GaAs and Al^Gai-^As. The re- 
sult is a quantum well of width 17.5 nm in the middle 
GaAs layer. Lateral confinement within that GaAs layer 
is due to the inhomogeneous electric potential generated 
by the top gate. In the original setup, a 30 nm thick 
GaAs cylindrical cap over the center of the dot created 
that inhomogeneity, while in our modified version, the 
different voltages on the inner and outer discs both play 
a role in the lateral confinement. The inhomogeneity is 
such that the potential energy of an electron V (= — e x $ 
where $ is the gate voltage) is lower near the axis of the 
dot so that electrons are attracted to the region under 
the inner disc or cap. 

We determine the potential profile within the dot by 
numerically solving the Poisson equation 



dl)V 





(_ e ) x -Sps 

v I Ke 



< z < 105 nm 
105 < z < 140 nm 



(39) 



applying the methods and parameters used in Ref. [33| to 
match experimental observed capacitance spectroscopy 
peaks of Ref. Here, z is measured from the top of 
the substrate, so that the region 105 nm < z < 140 nm 
is the charged blocking barrier. In the Poisson equation, 
K = 12.85 is the dielectric constant for undoped GaAs, 
and the charge density in the blocking barrier is taken to 
be no = 4.62 x I0 17 cm' 3 . We take the Fermi energy of 
the bottom electrode as the energy reference, and specify 
Vi n and V ou t at the top surface of the cylinder. Thus, 
at z = 140 nm, Vt n is taken to be a positive constant 
potential energy in the inner disk of radius Ri n and V ou t 
is a positive constant potential energy on the outer disc 
between Ri n and R ou t (to create electron confinement we 
must have Vi n < V ou t)- The potentials are defined to be 
those at the edge of the blocking barrier, z = 140 nm, 
so that Vin and V out include the Schottky potential of 
about 0.65 eV at the metal gate semiconductor interface, 
as well as the drop due to the GaAs cap under the inner 
disk. The boundary condition on the curved outer surface 
of the cylindrical dot is set by assuming a large value 
of R ut, so that the electric field is essentially parallel 
to the curved surface and aligned along the z-direction. 
We solve the ID Poisson equation in the z-direction to 
determine the potential profile on this outer surface at 
r = Rout 



V(R ou t,z) 



B lZ 

2 2 

ripe z 
' Ke 2 ~* 



-B 2 

Axz 



An 



< 105 nm 
> 105 nm 



(40) 



Applying the boundary conditions at the top (V — 
V ut) and bottom surface (V = 0) and matching 
solutions at the bottom edge of the blocking bar- 
rier layer at z = 105 nm determines the constants 
Ai = {Vout - 9.9684)/140 eV/nm,Bi = {V ou t - 
0.3987)/140 eV/nm, A 2 = 3.5886 eV, and B 2 = 0. 

With the value of the potential on the boundary deter- 
mined, the full 2D Poisson equation in r, z is solved using 
FISHPACK FORTRAN routines for the solution of sep- 
arable elliptic partial differential equations 3 -. Figure 0] 




Z [nanometers] 
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FIG. 4: Potential energy of an electron in the structure de- 
picted in Fig. E^b), without the 0.22 eV offsets due to the 
AlGaAs layers. 
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FIG. 5: Radial and axial profiles of the potential along lines 
passing through the center of the quantum dot obtained from 
the solution of the Poisson equation for the structure depicted 
in Fig. |3jb). (a) Axial (z) profile along the cylinder axes 
(R=0) including the 0.22 eV band offset (b) Radial profile 
at z=90 nm where the potential is the lowest; the black dots 
correspond to the solution of the Poisson equation, and the 
continuous line is the parabolic fit. 



plots the shape of the potential without the 0.22 eV con- 
duction band offset of the AlGaAs; we took Ri n = 120 
nm and R ou t — 400 nm. Fig.[5]Ja) shows the computed z- 
profile along the cylinder axis at R = nm now with the 
0.22 eV AlGaAs offset included; this shows that the re- 
gion where electrons are trapped resembles a finite square 
well along the z-axis. Fig. [5^b) shows the computed 
profile, in the middle of the finite square well [z = 90 
nm) , in the radial direction. The solid line is a parabolic 
fit, which shows clearly that in the radial direction the 
trapped electrons feel a harmonic oscillator potential. 
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VII. AN ELECTRON-PAIR IN THE MODIFIED 
QUANTUM DOT 

The problem of two electrons in a quantum dot has 
been analyzed in detail in the literature. As we see from 
Fig. [5j the potential in our quantum dot is well approx- 
imated by a finite square well in the z-direction and a 
harmonic oscillator in the radial direction. The effective 
Hamiltonian operator for two electrons is then 



H = — 



2m" 



[V 2 rl +d 2 zl +V 2 r2 + d 2 z2 ]+V sq . well 



(ri - r 2 ) 2 + (zi - z 2 ) 2 2 

In the numerical estimates that follow, we neglect 
the z— dependance in the Coulomb interaction in or- 
der to separate the radial and axial directions in the 
Schrodinger equation. Quantitatively, the neglect of this 
term in the denominator of the Coulomb interaction will 
cause a slight overestimate of the Coulomb interaction 
energy, an effect that can be diminished if the dot is 
significantly narrower along the z-axis than in the ra- 
dial direction. However, qualitatively this approximation 
makes no difference for the demonstration of quantum 
pumping feasibility. 

The total energy of an electron in the dot can now be 
separated into two parts E = E z +Er, an axial quantum 
well energy Ez which is well approximated by a finite 
square well and a transverse energy Er that can be well 
approximated using a harmonic oscillator. 



A. Axial energy, E z 

The size of the quantum well is L — 17.5 nm. and the 
height is given by U = \U \ + 0.22 cV, where 0.22eV is 
the band offset between GaAs and Al^Gai-^As and Uq 
is the depth of the well below the Fermi energy of the 
bottom lead, our reference energy. For an infinite well of 



the same width we get: E^° = 



2h z 

in* L 2 



(§) = 18.31meV 



taking m* = 0.067m for undoped GaAs. For a finite 
square well, the ground state energy is determined by 



2ftV 
m*L 2 



(7.43meV) X X \ X tan(x) = VP* 



with P 2 = 134.8 x C/(in eV) where U is the well depth. 
We found that as U varies over the range (0.22, 0.26) 
eV, E z remains in the range (13.0, 13.4)meV. Since it 
turns out that Uq changes by less than O.OleV when we 
vary the parameters V ou t and Vi n , the value of Ez , which 
is defined relative to the bottom of the well, changes by 
less than ~ 1%. Therefore, in our calculations, the vari- 
ation of Ez relative to the bottom of the well can be 
neglected. As Fig. EJa) shows, the bottom of the well is 
not completely flat, varying from 0.22eV to \U \ +0.22eV. 
However, for the same reason, this slope has little effect 
on the constancy of Ez- 



B. Radial energy, Er 

The transverse or radial Hamiltonian can be separated 
into center-of-mass (R) and relative coordinates (r)2£ 



H, 



- 2 D 2 

-7 R 



-2V z r + - 7 V + 



(42) 



where energies are scaled in Rydbergs Rd = 2m * a 'i 
and the length scale is ajj — h 2 (AireoK)/(m*e 2 ) and 
7 = 2(a_o/Zo) 2 where lo = \Jh/m*LO is the oscillator 
length. In these units, the harmonic oscillator energy 
spacing, Tito, is given by RdJ- Excluding the Coulomb 
interaction the solutions are those for a 2-dimensional 
harmonic oscillator with circular symmetry with energy 

E R - E mt = (2N + 1 + |Af |)7 + (2n + 1 + |m|) 7 . (43) 

The parity of the spatial wavefunction is determined by 
the relative coordinate quantum number m. For m = 0, 
the spatial wavefunction is even under particle change 
(r — > — r), so the state must have an antisymmetric 
spin wavefunction, a singlet. For m — 1 the spatial 
wavefunction is odd under particle interchange, so the 
state must have a symmetric spin wavefunction, a triplet. 
We conclude that the lowest energy states are the sin- 
glet with (N,M,n,m) = (0,0,0,0) and the triplet with 
(N,M,n,m) = (0,0,0,1). These have energies 2huj and 
37iw respectively without including the Coulomb interac- 
tion energy E int . 

The interaction energy Ei nt has been calculated for 2D 
circular quantum wells in the literatur e) 35 ' 36 for values of 
7 ~ 0.1 — 1. We use those results for E int to determine 
the total radial energy Er. 



VIII. EXPERIMENTAL FEASIBILITY OF 
PUMPING ENTANGLED ELECTRONS 

A. Desired scenario for singlet pumping 

We would like to vary the top gate potential energies 
Vi n and V ut m such a way that the shape of the lat- 
eral confining potential changes with the following con- 
straints: 

(i) At all times, the singlet state is energetically ac- 
cessible, but the triplet state is too high in energy to be 
occupied. 

(ii) The energy of a single electron in the dot with re- 
spect to the Fermi energy Ep = of the bottom electrode 
stays constant and far off resonance with Ep. This sup- 
presses single electron tunneling, which depends strongly 
on including the resonant region in the pumping cycled. 

(iii) There is significant variation of the Coulomb 
interaction between the two electrons: the singlet energy 
can be varied appreciably to permit pumping of pairs. 
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We present numerically computed results for two sam- 
ple configurations of the dot corresponding to different 
specific values of the top gate potential energies Vi n and 
V ou t which satisfy the above criteria. These two configu- 
rations could serve as points within a cyclic variation of 
the top gate potential energies. 

The lowest energy state spin singlet (with quantum 
numbers (N, M, n, m) = (0, 0, 0, 0)) and spin triplet (with 
quantum numbers (N, M,n,m) = (0,0,0,1)) have total 
energies: 

Is : £ (S) (w) = 2U + 2E Z + 2huj + E^} (w) 

2p : E (T ^ (w) = 2U + 2E Z + ihu + E$ (w) (44) 

In order to satisfy condition (i) we look for values of lu 



such that the total triplet energy lies above the Fermi 
level i.e. the reference energy in our calculations, while 
the E^ lies below: 

£ (S) H<0 £ (T) H>0 (45) 

In order to satisfy condition (ii) to maintain the sin- 
gle particle energy levels fixed with respect to the Fermi 
level, we have to simultaneously ensure that in any two 
configurations (1 and 2), we have: 

hwx + +E z = hu)2 + uP +E Z , (46) 
=>■ TvjJx - Tioj 2 = UP - Ui 1} . 



TABLE I: Summary of numerically computed parameters for two specific configurations of the potential profile of the proposed 
quantum dot, that would satisfy the requirements of singlet pumping. The two configurations would be part of a continuous 
and cyclic variation of the potential profile of the dot. 



Configuration 1 


Configuration 2 


Vout 
Vin 


1158.1 meV 
263.1 meV 


950.6 meV 
301.2 meV 


Single Electron level 


la = \Z%/m*u 
hu 
U 
Ez 

Uo + E z + hoj 


14.35 ran 
5.5 meV 
-22.9 meV 
13.2 meV 
-4.2 meV 


15.55 nm 
i.7meV 
-22.1 meV 
13.2 meV 
-4.2 meV 


Singlet interaction energy 
Triplet interaction energy 


E (S) 

P (T) 

E {T) - E (S) 


7.3 meV 
4.6 meV 
-2.7 meV 


6.6 meV 
4.2meV 
-2.4 meV 


Radial Singlet energy 
Radial Triplet energy 


E (S) 
^R 

^R 

^R 


18.4 meV 
21.2 meV 
2.8 meV 


16.0meV 
18.3 meV 
2.3 meV 


Total Singlet energy 
Total Triplet energy 


2U + 2E Z + E^ } 
2U + 2E Z + e£ } 


-0.5 meV 
+0.9 rneV 


-0.9 meV 
+0.25 meV 



B. Experimental parameters for singlet pumping 

We take the dimensions of our quantum dot to be simi- 
lar to those of Ref. [32|, except that our dot has a slightly 
smaller radius of 120 nm for the inner top gate while 
estimates^ found that in the experiment the radius of 
the GaAs cap was about 200 nm. 

We explicitly identify two sample locations in param- 
eter space that satisfy all of our conditions. (To find 
a complete pumping path through parameter space is 
straightforward since one has two parameters V ou t and 
Vi n and only one precise quantitative constraint (ii) 
above). Say the two sample locations have radial oscilla- 



tor energies Too = 5.5 meV and Tilj — 4.7 eV (correspond- 
ing to 7 = 1 and 7 = 0.85 respectively). Then the well 
depths must compensate such that the single electron en- 
ergy (|46p remains constant. For example, the locations 
could satisfy Uq = —22.9 meV and Uq = —22.1 meV so 
that the single electron energy is held at —4.2 meV, con- 
stant and strongly off resonance with the Fermi energy. 
At these values of hu> and Uq, the quantum dot geom- 
etry is such that the interaction energies^ for singlets 
are 7.3 meV and 6.6 meV respectively and for triplets 
4.6 meV and 4.2 meV, respectively. Then equation (|44p 
implies that the total energies for the singlets are —0.5 
meV and —0.9 meV for the two locations in parameter 
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space while the total energies of the triplets are 0.9 meV 
and 0.25 meV respectively. Thus, the singlet states are 
energetically accessible while the triplet states are not. 

Our simulation shows that the two configura- 
tions are realized in the dot at (V ou t,Vin) = 
(1158.1meV,263.1meV) and (V out ,V in ) 
(950.6meV, 301.2meV), from our numerical compu- 
tation, which can be rounded off to the appropriate 
significant figures. Figure^a) shows the radial profile in 
the dot for the two configurations. Configuration 1 has 
a higher difference between Vi n and V ou t, leading to a 
larger value of us and causing tighter lateral confinement. 
Configuration 2 has a lower difference between V% n and 
Vout leading to smaller value of us. However, because Vi n 
is higher than in configuration 1, the well is not as deep. 
This can be clearly seen in Fig. [6]Jb) where the profiles 
for the two configurations cross due to the compensating 
effects of increased us and lowered Uo- 

Table I summarizes the results of our simulation. Ener- 
getically the singlets are accessible in both configurations, 
while triplets are inaccessible (point (i) above). The en- 
ergy of the single electron level in the dot is fixed at 
— 4.2meV significantly off-resonant with Ep = (point 

(ii) above). The variation of the Coulomb interaction is 
about 11% between the two singlet configurations (point 

(iii) above). 



IX. CONCLUSIONS 

In this paper we address several issues related to de- 
veloping a quantum singlet pump. We will conclude by 
providing a summary of our main results. We have con- 
sidered the problem of applying the mechanism of adi- 
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FIG. 6: (a) Radial profiles through the lowest potential region 
of the dot for two distinct set of top gate potential energies: 
continuous line — > (V ou t, Vin) = (1158.1 meV, 263.1 meV) and 
dotted line -> (V ou t,V in ) = (950.6 meV, 301.2 meV). They 
correspond to the two configurations described in Table I. 
(b) Expanded view of the lowest energy region. The shift 
Uo between the two configurations is Tiusi — T1W2, keeping the 
single particle ground state energy unchanged. 



abatic quantum pumping to generate a flow of singlet 
pairs of electrons, while suppressing the flow of triplets 
and uncorrelated single particles. The first challenge was 
to find the appropriate theoretical description. We first 
identified an appropriate definition of the current for sin- 
glet pairs in analogy with the current of a stream of un- 
correlated electrons, by using the reduced two-particle 
density matrix. We then showed how in the presence of 
two-body interactions, the evolution of the many elec- 
tron system can be effectively described by the evolution 
of a two-particle state particularly when the interaction 
is spatially localized. We confirmed that our definition of 
current gives a non-zero current in the absence of bias or 
time variation. We then wrote an adiabatic perturbation 
expansion for the time varying two particle states, where 
the rate of change of the interaction dV/dt is assumed 
slow compared to the time scale of the dwell time of the 
particles in the interaction region. Using the assump- 
tion of adiabaticity, we computed the pumped current to 
first order in the time derivative (or equivalently the fre- 
quency of the time- varying interaction) . The interaction 
is assumed to affect only singlets and is non-vanishing 
only in a finite region, but otherwise completely general. 
By using several identities and relations for two particle 
Green's function, that we present in the appendices, we 
reduced the expression for the current to a compact and 
relatively simple form involving only the instantaneous 
two particle Green's functions. The current is seen to 
have an additional transient term which we show van- 
ishes identically for a complete cycle, the term being a 
total time derivative. 

Having established a general but simple expression for 
the pumped singlet current due to the action of a time- 
varying local two body interaction, we apply it to a spe- 
cific model. We consider a lattice, where the interaction 
acts and varies in time only at two of the sites; the in- 
teraction acts only when both particles are localized on 
the same site, while the two sites are sufficiently sepa- 
rated that interaction among an electron at one site with 
one at the other can be neglected. If the sites are well 
localized in space, the effect of the interaction on triplets 
is automatically suppressed due to vanishing or at least 
substantial suppression of the triplet spatial wavefunc- 
tion at each of the two sites due to the Pauli principle. 
We computed the singlet current for the model as an ex- 
act analytical expression, which has been confirmed to 
have the same form for discrete and continuous models. 

Finally we show that this model can be implemented 
using a chain of quantum dots. We take the design spec- 
ifications of a quantum dot that had been fabricated and 
studied in the lab and propose one minor modification 
that would introduce two concentric top gates with in- 
dependently controllable voltages, something that can be 
achieved without much difficulty with the methods of fab- 
rication available currently. We computed the potential 
profile of such a dot with available experimental param- 
eters and computed the energy of two electrons in such 
a dot. We showed that by varying the two gate volt- 
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ages, (i) significant singlet current can be generated (ii) 
triplets can be made energetically inaccessible for states 
in the dot, and therefore affected much less by the varia- 
tion of the interaction (iii) the single particle states can 
be maintained at the same energy and far off-resonance, 
thereby suppressing current of uncorrclatcd single elec- 
trons. Thus in effect such a dot can be used in a chain 
to implement our model for generation of a measurable 
singlet current. 
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APPENDIX A: GREENS FUNCTIONS: 
REPRESENTATIONS AND IDENTITIES 

The Green's function for the potential-free time- 
independent Schrodinger equation is defined by 



[E+^-V 2 }G (x',x";E)=S(x'-x") 



(Al) 



where x = {x\,x 2 , ••• ,x n } denotes the spatial degrees 
of freedom, the total energy E shared among them is as- 
sumed real, and the subscript in Go indicates the absence 
of any potential. The degrees of freedom have equiva- 
lent interpretations as n spatial dimensions or as coor- 
dinates of n individual particles in one dimension. We 
study interacting particles in ID, but we use 'n-particle' 
or 'nD' interchangeably in referring to the Green's func- 
tions. The retarded and advanced Green's functions will 
be denoted using superscripts G + and G~. When the 
Green's functions obey G$(x', x"; E) = [G (x' , x"; E)}* , 
we leave out the superscripts, G + = G and G~ = G*. 

Our interest being in two body interactions we only 
need single particle (or ID) and two particle (or 2D) 
Green's functions. The uppercase G(x' , x"; E) will be re- 
served for 2D Green's function with x = {xi, x 2 }, and ID 
Green's functions will be distinguished by using lowercase 
go{x,x';E). Single-particle eigenstates are likewise de- 
noted by lowercase 4>k{x); expansion in such eigenstates 
readily establishes the following useful identities: 



£ 
£ 



dx<t>l(x)g (x,x';E) 

g (x',x";E-e) 

e- Ep + ir] 
g (x',x";E-e) 



<t>%{x') 



de 

> 

de 



E - E k + ir/ 
-2iirg (x',x";E- E p ) 



(A2) 



En 



IT] 



= 



dxgl(x,x'\E 1 )g (x,x";E 2 ) = 

gZ(x",x';E 1 )-g {x',x";E 2 ) 



Ei - Ei + ii] 



/oo 
dxg (x,x , ]E 1 )g (x,x":E 2 ) = 
-00 



g (x" ,x';Ei) - g (x' , x" ; E 2 ) e 2 ^e 1 d , „ 
p p ~* ~7T~ 9( - x ' x ' El > 

£j 2 — &\ OE! 

For scattering problems the appropriate eigenstates are 
plane waves which in ID are <pk {x) = e lkx ; they provide 
a coordinate representation for the free ID Green's func- 
tion for real energies E 



9o(* 



',x";E)= f 

J —I 



dk 



a ik{x' — x") 



2irE- {h 2 k 2 /2m) ± ir] 



(A3) 



= T 



i6(E)e 



±1 



j2£E\x'-x"\ 



g i -E)e-^ F ^'-' 



h 2 



^2mE/h 2 



■sj2m\E\/h 2 



The difference of the retarded and the advanced Green's 
functions gives the two-point correlation function 

p(x',x";E) = '' 



-2ni 



[g+(x',x";E)-g (x",x';E)](A4) 



dk 



6(E) I —S(E~(h 2 k 2 /2m))e lk{x - x '"> 
J-00 2?r 

Performing the momentum integral demonstrates consis- 
tency with the definition of gf . In the case of equal 
coordinate arguments a different contour integral is in- 
volved, but the end result agrees with simply setting 
x 1 = x" in the above expressions; p(x' , x'; E) defines the 
one-dimensional density of states. 

The 2D free Green's function can be written in terms 
of the ID Green's functions 



G$(x'x";E) 



(A5) 



±i 
2^ 



/CO 
dego {x[ , x'{ ; E - e)g^ (x' 2 , x 2 ; e) 
-00 



Using an eigenfunction expansion in Cartesian coordi- 
nates and integrating out one momentum component 
gives 

G±(x'x";E)= (A6) 
'■\fW dk2cos[k{x' 2 - 40]e ±i \ /2ro - E/ft2 - fe2|x ' 1 - x " 1 



Jo 



2tt 



2mE/h 2 - k 2 



dk 2cos[fc(4 - x »)] e -V k2 - 2mE / h2 \ x 'i- x i\ 



2mE 2lT 

K 2 



-^k 2 ~2mE/Ti 2 



The first term has both real and imaginary parts and is 
always regular (R) as a function of the coordinate ar- 
guments, while second term is always real and is sin- 
gular (S) when x' — x"; so we name the two terms 
Gq R (x'x";E) and G a t s (x'x";E) respectively. 
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APPENDIX B: INTEGRALS USED IN 
COMPUTING CURRENT 

We will apply the results of the preceding appendix 
to evaluate the generic expressions involving two particle 
or 2D Green's functions required in computing current. 
The 2D Green's functions below are at energy E, not 
explicitly shown for the sake of compact notation. 



1. Integral with conjugate Green's functions 

We encounter expressions involving a pair of conjugate 
2D Green's functions when computing the current: 



f°° d 
I g *g(x',x")= J dx 2 G* (x,x') — Go{x,x") (Bl) 

By expressing the 2D Green's functions above in terms 
of ID Green's functions as in Eq. (|A5|) . and then using 
the identities in Eq. (1A2|) we can reduce it to the form 



I g *g(xi,x',x") = — / de[g{x' 2 ,x 2 ;e) - g* {x' 2 , x 2 ; e)] 
xg* ( Xl ,x[;E-e) ^-g(x x , x'{; E - e). (B2) 

OXi 

Inserting the expressions for ID Green's functions from 
Eqs. (|A3| and (|A4|) leads to an explicit functional form 
that is piecewise continuous, in which the free variable 
xi determines the boundaries of continuity; thus for the 
exterior region |a;i| > \x'i\, \x'{\, which is relevant to us, 
the expression is 

I G * G (x 1 ,x',x")= (B3) 
dk 2cos[k{x 2 — x 2 )\e v R 



± 



2tt 



Kj2mE/h 2 -k 2 



j; o r7 I I IIW - \/ k2 -TT E \' 2x i-( x 'l+ x l)\ 

dk 2cos[k(x 2 ~ x 2 )\e V ^ 



2tt 



k 2 - 2mE/h 2 



where the + sign applies for x\ > x'i,x'{ and the — sign 
for x\ < x'^x'l- We can exploit the similarity of this 
expression with G^(x'x") in Eq. (|A6[) because we always 
encounter it in the functional combination 



F x I G » G (xi ) x' ) x") 
Im I F 



d 

dx 2 Gl{x, x')- — Go(x, x' ) 
i ox i 



f°° d 
+F* / dx 2 G* {x,x") — Go{x,x') 

J-oo OX\ 



where F is a complex-valued function and x' <-> x" are 
exchanged between the two terms; the second term in 



Eqs. (|B3p . being (i) manifestly real and (ii) unaffected by 
the exchange x! «-> x", would not contribute to Eq. (|B4|) . 
Therefore, in this particular combination the integrals 
can be replaced by G^ R (x'x"). We need to note that the 
exponential in Eq. (|A6|) contains the absolute difference 
of the coordinate arguments while Eq. (|B3[) does not; 



this in effect determines the choice of G 0R depending on 
whether x' > x" or x' < x" : 

For xi > x' 1>X '{, F x I g >g(x', x") = (B5) 

glm {-[FG 0R (x', x") - F*G* 0R (x', - x'{) 

n 

+[FG* 0R (x',x") - F*Gor(x',x")]6(x'; ~ x[)}} 
For xi < x'i,Xi, F x Ig' G (x, X ) = 

■^\m{~[FGl R {x ,x") - F*Gqr(x ,x )]6(x t - x x )] 

+ [FGor{x',x") - F*G* 0R (x',x")]6{x'{ - x[)]} 
It is also useful to note that the imaginary part satisfies 

WW*',*")} = {^frf^'l Xl> /> ,X l 

[ +^Im{G (a;',a;")} xi <x' 1 ,x'{ 

which helps simplify the CclSG X — X . when the first 
terms in both Eqs. (|B3j) and (|A6j) become imaginary 

F x I g *g(xi,x',x') 

= f -^lm{G (x',x')}Re(F) x x > x\ 
\+§lm{G (x',x')}Re(F) x x <x' x ' 

2. Integral with like Green's functions 

We also encounter expressions with a pair of similar 
2D Green's functions 

I GG {xi,x',x")= / dx 2 G (x,x')—G {x,x") (B7) 
J-oo oxi 

which, by using Eq. (|A5|) and the identities in Eq. (|A2[) . 
can be reduced to the form 
i f°° 

= — J deg(x' 2 ,x 2 ;e)g(xi,x' 1 ;E - e) 

x^-g( Xl ,x'{;E-e) (B8) 
ax i 

Then the expressions for the ID Green's functions in 
Eq. (|A3[) provide explicit functional forms which, as in 
the previous case, depends on the value of x\. For the 
exterior region \xi \ > \x'i\, \x'{\ this is 



(B4) I GG (xi,x' ,x") = ± 



2- 



deg e (x 2 ,x 2 ) 



(B9) 



Jy / 2m{E-e)/ti 2 \2x 1 -(x' 1 +x")\ 



i^\l2m{E-e)/tf 



-8(E-e) 



-y / 2m(e-E)/h 2 \2x 1 -(x' 1 +x")\ 



h 1 



^2m(e-E)/h 2 



(e-E) 



15 



where the + sign applies for x\ > x[,x'{ and the — sign 
for x\ < x\ , x'{ . This integral is also of interest for inte- 
rior region when |a;i| is between \x'i\ and \x'{\: 



t ( -i -n\ I +¥?Go(x' , x") x'i > x\ > x'{, , 
I GG {xi,x ,x ) = < n ) ' 1 J(B10) 
I -§iGo(x\x") x'{ > xi > x\ 



The following combination is relevant 
F x I GG (xi,x',x") = 

ImlFj dx 2 G {x,x')—G (x,x") 



(Bll) 



+F* / dx 2 G* {x,x")—G* (x,x') 

, OXi 



For the exterior region I GG (x±, x', x") is invariant un- 
der exchange of coordinate arguments, so the expres- 
sion above vanishes, F x I GG (xi,x' ,x") = 0, for the 
exterior region \x\\ > \x / 1 \,\x / [\. In the interior region 
Igg{x\, x', x") simply changes sign under such an ex- 
change, so that for 



APPENDIX C: AN EFFECTIVE TWO PARTICLE 
STATE FOR THE MANY BODY SYSTEM 



As soon as we consider pair interaction, in principle, we 
have to allow for interaction between all possible pairs in 
the system. In this Appendix we reduce a full many body 
system with pair interaction to an effective description in 
terms of a pair of interacting particles both moving in the 
background of the interaction field due to all the other 
electrons in the system. This allows us to describe the 
pumping of singlets in terms of the evolution of a singlet 
wavefunction very much like the pumping of individual 
electrons by a single particle state. 

Using the composite notation for position and intro- 
duced in Sec. II, X = (x, a), the Heisenberg equation of 
motion of a Fermion creation operator with the Hamil- 
tonian ([7|) is 

ihd t ^(X,t) = h(x)$(X,t) 

dX'ft(X, i)ft(X', t)V(x, x')tp(X', t). (CI) 



F x I GG (x',x") = ±^lm{FG (x',x") - c.c.} 

= f +2$Im{FG (x',x")} , x[ > x x > x'{ 
\ -2§lm{FG (x',x")} , x'{ > Xl > x[ 



(B12) 



Consider (E , N\fo (Xx,t)fo (X 2 ,t) \E 2 ,N- 2), a matrix 
element where \E ,N) is the N particle ground state of 
the system, and \E 2 ,N — 2) is an energy eigenstate with 
N— 2 particles. Using Eq. (|C1[) . we find that the equation 
of motion of this matrix element is 



ind t (E ,N\^(X 1 ,t)^(X 2 ,t)\E 2 ,N - 2) = {h{ Xl ) + h(x 2 ))(E Q , N\ft(X u t)ft(X 2 , t)\E 2 , N - 2) 

+ J dX'V(x 1 ,x')(E ,N\^(X 1 ,t)^(X',t)^(X',t)^(X 2 ,t)\E 2 ,N-2) 

+ f dX'V(x 2 ,x')(E ,N\ft{X u t)ft(X 2 , t)ft{X', t)iJ;(X', t)\E 2 ,N - 2) (C2) 



Since both \Eo, N) and \E 2 ,N—2) are energy eigenstates, 
we can replace ihdt — > E = E 2 — Eq on the left hand side 
of the above equations. Furthermore, 

fi(X',t)j>(X\t)j>\X 3 ,t) (C3) 
= S(X' - X 2 )$ (X', t) + ft(X 2 ,t)ft (X', t)${X\ t) 

due to the Fermion commutation relations; here we 
denoted S(X' — X 2 ) = S a ^ a2 S(r' — r 2 ). Making a 
mean-field approximation, with the density of back- 
ground electrons at X' defined by, n(X',t) = (E ,N — 
2\^ (X' ,t)ijj(X' ,t)\E Q , N - 2) , we define the two particle 
state 

* E (X 1 ,X 2 ,t) = (E Q , N\ip\Xi , t) i>\X 2 , t) | E 2 , N - 2) (CM) 
and bring (|C2j) into the form 



E^ E {X l7 X 2 ,t) (C5) 
= [h(xi) + h(x 2 ) + V(x 1 ,x 2 )]^ E (Xi,X 2 ,t) 

+ J dX'[V{x 1 ,x') + V(x 2l x')}n(X\t)^ E {X 1 ,X 2 ,t) 

The last line gives the mean-field influence of the 
rest of the electrons with the particular pair of elec- 
trons at Xi and X 2 . We can include this in the one 
body potential that each particle experiences: h[x) + 
/ dX'V(x,x')n(X' ,t) — > h(x). Since the two-body inter- 
action V does not affect spin, we can factorize out the 
spin part of the wave function. Considering specifically 
the singlet subspace and by introducing the parametriza- 
tion in terms of the single particle momentum labels 
k% , k 2 we are led to Eq (jHJl . 
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